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The background field method for measuring the electric polarizabiUty of the neutron is adapted to 
the dynamical quark case, resulting in the calculation of (certain space-time integrals over) three- 
and four-point functions. Particular care is taken to disentangle polarizability effects from the 
effects of subjecting the neutron to a constant background gauge field; such a field is not a pure 
gauge on a finite lattice and engenders a mass shift of its own. At a pion mass of m^j; = 759 MeV, 
a small, slightly negative electric polarizability is found for the neutron. 
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1. Introduction 

The response of hadrons to external electromagnetic fields is characterized by their polar- 
izabilities. Polarizabilities manifest themselves, e.g., in soft Compton scattering; the incoming 
electromagnetic wave distorts the hadron, and this distortion in turn leads to a modified scattered 
wave. Accordingly, polarizabilities appear as effects of second order in the external field in the 
effective hadron Hamiltonian. In a low-energy expansion of that Hamiltonian, the leading term 
embodies the hadron mass shift in constant electric and magnetic fields E and B, 



with the static electric and magnetic polarizabilities a and j8, respectively. Carrying the low-energy 
expansion to subsequent orders, dependences on particle spin and electromagnetic field gradients 
appear [|l]]. 

The work presented here is concerned specifically with the static electric polarizability a of the 
neutron. There are two chief new elements to this investigation. On the one hand, dynamical quark 
ensembles are used; on the other hand, the role of constant background gauge fields is recognized 
and their effects are properly taken into account. As far as the former issue is concerned, it should 
be noted that, when introducing an explicit external electromagnetic field to compute polarizabili- 
ties on the lattice, the transition from quenched to dynamical quark ensembles entails more than just 
the usual increased computational cost of generating the ensemble. In the quenched approxima- 
tion, the backreaction of the quark fields on the gluonic degrees of freedom is truncated; thus, the 
external electromagnetic field, which in turn only couples to the quark fields, does not influence the 
gauge ensemble. As a result, one can generate the ensemble in the absence of the background field 
and introduce the latter a posteriori by a suitable modification of the link variables. By contrast, 
in the dynamical quark case, at first sight, one would need to recompute a new gauge ensemble 
for every specific external field considered, a prohibitively expensive approach if directly applied. 
Previous quenched investigations |^, ||, ^ were thus able to take advantage of considerable sim- 
plifications not available for dynamical quarks. The dynamical quark case necessitates significant 
modifications of the computational scheme in order to become tractable, cf . section |2| below. 

Concerning the issue of constant background gauge fields, consider introducing a constant 
electric field in the 3-direction via the gauge field 



For given electric field E, one still has a continuum of choices for Iq or, equivalently, A. In a setting 
in which the spatial directions are infinite, this is merely a gauge freedom; A can be arbitrarily 
shifted by gauge transformations. However, if the spatial directions are finite, such as in a lattice 
calculation, this gauge freedom is reduced by the need to preserve the boundary conditions; as 
a consequence, physics in general depends on the value of A. As a simple example, consider a 
charged particle on a circle of length L with periodic boundary conditions in a constant gauge field 
A. The Hamiltonian is H = (p — A)^ and the energy levels are = {2nn/L — A)^, where n is any 
whole number. In particular, the ground state energy as a function of A is plotted in Fig. Only a 
discrete symmetry under shifts A — > A + 27r/L remains; A can be interpreted as a Bloch momentum 
of the charged particle. 




(1.1) 



Asit) =E{t-to) =A + Et . 



(1.2) 
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Figure 1 : Ground state energy of a charged particle on a circle of length L with periodic boundary conditions, 
as a function of the constant gauge field A to which it is subjected, cf. main text. 



Similarly, the mass shift of a neutron in the external field ( |1.2[ ) is a function of both the pa- 
rameters E and A. The dependence on A is a dominant effect which must disentangled from the 
polarizability in presence of the electric field E. 



2. Perturbative expansion 

As mentioned above, in the presence of dynamical quarks, the gauge ensemble depends on 
the specific external electromagnetic field under consideration. However, it would be prohibitively 
expensive to recompute the ensemble for each background field. Relief is provided by the fact that 
it is not necessary to be able to evaluate the neutron mass shift in arbitrary backgrounds; after all, 
the electric polarizability is given specifically by the second Taylor coefficient of that mass shift as 
a function of the electric field, cf. (|lT|). To evaluate a Taylor coefficient, it is sufficient to consider 
only infinitesimal external fields, and in this case, one can correspondingly expand the neutron 



two-point function determining the mass shift. Introducing the external field A3, cf. (|1.2|), into the 
link variables as an additional U{1) phase, 

U3 — > exp (^i J dx3 • A3^ .U3 = {1 + iaA3 - a^Aj/l + ...)• t/a , (2.1) 



one obtains an action 



S — So + Sext 



(2.2) 



composed of the QCD action in the absence of the external field, ^o, and the coupling to the external 
field 



Zv^I, Yix) [(/aA3 - a^Al/2) ■ U^ix) • (-1 + 73) • w{^ + e^)) 

+{iaA3 + a^Al/1) -Ulix- €3) ■ {\ + Y3) ■ wi^- e^)) 



(2.3) 



which has been written here specifically for Wilson-type four-dimensional quark fields. The par- 
ticular computational scheme adopted in the following uses domain wall fermions in the valence 
sector; however, quark propagators are evaluated starting from sources which only have support on 
the four-dimensional domain walls, and they are projected back onto four dimensions at the sink. 
This is a compromise due to hardware constraints on the storage of propagators; correspondingly, 
the coupling to the external field is effected via the four-dimensional projected fields i/A, as written 
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in Concomitantly, it is necessary to allow for the renormalization factor zv on quark bilinears 
to compensate for the effect of the projection. This factor was determined in [§] by evaluating the 



number of valence quarks in the neutron, yielding zv = 1.12 ±0.1 2. Inserting the action ( [2.2| ) into 
the neutron two-point function, one has 

{Nfi{y)Na{x)) = J[DU][DY]m]exp{-So) {l-S,^ + Sij2 + ...)Np{y)Na{x) / 

J [DU] [Dy] [Dv/] exp(-5o) (l - 5,,, + S^,^/2 + ...) (2.4) 

To obtain specifically the part of the two-point function which is quadratic in the external field, one 
therefore must calculate neutron two-point functions with either two additional insertions of the 
vertex linear in the external field, cf. (|2.3[), or one additional insertion of the vertex quadratic in the 



external field, cf. (2.3), as displayed in Fig. 



O 

n n n n n n n 

Figure 2: Diagrams of second order in the external field resulting from inserting the electromagnetic vertices 
generated by the interaction term Sext into a neutron two-point function. Crosses denote vertices linear in the 
background field, circled crosses vertices quadratic in the external field. 



The space-time positions of the insertions are to be integrated over according to ( |2.3[ ); thus, one 
is effectively calculating certain space-time integrals over three- and four-point functions. These 
functions are evaluated in the unperturbed ensemble generated by the action Sq. The two right-most 
diagrams involving the vertex quadratic in the external field act as contact terms renormalizing the 



second and fourth diagrams from the left, respectively. As usual, the denominator of (2.4) effects 
the subtraction of statistically disconnected pieces from diagrams containing disconnected quark 
loops. 

3. Extracting the mass shift 

Focusing on the zero-momentum, unpolarized component of the neutron two-point function, 

G{p = 0,t)='L^Tr(^^^{Niy)Nix))^ W expi-mt) (3.1) 

one observes an exponential decay in time at a rate governed by the neutron mass m; W charac- 
terizes the overlap between neutron sources/sinks and the neutron ground state wave function. In 
the setting discussed here, some care must be exercised in interpreting the two-point function [S]. 



The Hamiltonian in the presence of the external field ( |1.2[ ) is time-dependent, since a translation 
in time is equivalent to a shift in the constant gauge field A, cf. ([T^); on the other hand, a change 
in A modifies the physical spectrum. For infinitesimal external fields, this time dependence can be 
treated adiabatically, i.e., the exponential decay of the two-point function can still be interpreted in 
terms of the neutron mass, which, however, will continuously adjust as the external field changes. 
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Expanding the overlap W and the neutron mass m in terms of the external field parameters A and 
E, 



W = Wo + w(^)(A,£') + WP)(A,£) + , 
m = mQ + nP'\A,E) + . . . 



(3.2) 
(3.3) 



where the superscripts denote the order in the external parameters, the second-order part of 
G{p = 0,t) behaves as 



\ Wo 



(3.4) 



Note that, in (3.3), use has been made of the fact that there can be no first-order term in the neutron 



mass in the absence of a nontrivial 0-angle. Thus, dividing (3.4) by the unperturbed correlator 
Wo exp(— mo?) yields a ratio whose slope as a function of time, S{A,E), corresponds, up to a 
minus sign, to the mass shift of the neutron in the presence of the external electromagnetic field. 



-S{A,E) 



m 



(2) 



(A, £) = m^A^ + n4^AE + m^'^E^ 



(3.5) 



Here, again, care must be taken in interpreting the results. As mentioned above, the spectrum of 
the system adjusts adiabatically to the external field, i.e., m^^' and W^"^^ themselves in general are 
time-dependent. This complicates the interpretation of S in terms of the mass shift ni^^^. However, 
at the particular point where the physical spectrum is stationary in A, the time dependence of m^^^ 
and W^^^ is of at least quadratic order. This is due to the equivalence between shifts in A and shifts 



in time, cf. (1.2); stationarity in A implies stationarity in time. Thus, at the stationary point only, 



which in practice is identified by seeking out the extremum of the slope S{A,E) as a function of A, 
one can indeed directly interpret that slope as giving = —S. 

In turn, as discussed in more detail in [^, the part of the neutron mass shift due specifically 
to its electric polarizability can be isolated from the combined dependence of m^^' on E and A by 
seeking out the extremum of m^'^\A,E) (or, equivalently, S{A,E)) as a function of A. The mass 
shift obtained at the extremum then yields the electric polarizability via 



a 



(3.6) 



cf. (1.1). To carry out this program in practice, one thus needs measurements at three different 



values of A (or, equivalently, to) in (|1.2D, while E is fixed. This then defines a parabola in to, 
permitting the extraction of the extremum as a function of to- 



4. Measurements 



The program described above was put into practice using 99 configurations from one of the dy- 
namical 2-1-1 flavor Asqtad quark ensembles provided by the MILC collaboration [^]. Specifically, 
the SU (3) flavor-symmetric ensemble with quark masses am„ = amj = arris = 0.05 was utilized, 
where a = 0.124fm denotes the lattice spacing. This choice of quark masses corresponds to a pion 
mass of niji = 759 MeV. On the other hand, as already mentioned in section^ in the valence sector. 
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connected 
diagrams 


all 
diagrams 


S/(a^E^) 
?Q = - IQa 


0.46(18) 


0.26(26) 


S/{a^E^) 
to = 


0.000(16) 


-0.033(43) 


?o = da 


-0.017(3) 


-0.037(25) 


-m(2)/(a3£2) 
(extremum) 


-0.034(6) 


-0.052(24) 



Table 1: Measurements of the slope S{A,E), together with a determination of its extremum as a function of 
A. The electric field E is cast in Gaussian units. Errors were obtained by the jackknife method, including 
the fluctuations of the renormalization constant zv- 



domain wall quarks were used; this hybrid scheme has been previously employed in ||], where 
also details concerning parameter tuning in this approach can be found. 

Table |l] displays the measurements obtained for the slope S{A,E), together with a determina- 
tion of its extremum as a function of A, which yields (up to a minus sign) the mass shift due to the 
electric polarizability of the neutron. In all measurements, the neutron source was placed at f = 
and the slope S{A,E) was extracted by fitting in the range Aa <t < lOa. The disconnected traces 
were evaluated using stochastic estimation with 120 complex Z(2) stochastic sources (240 for the 



case Iq = —10a, which displayed particularly strong statistical fluctuations). Using ( |3.6| ) to obtain 



the electric polarizability of the neutron from the last line of Table |], one finally arrives at 

a = (-2.0±0.9)-10"'*fm^ (4.1) 

at a pion mass of niji = 759 MeV. 



5. Discussion 

Compared to the experimental value reported by the Particle Data Group [^, namely, a 



(1 1.6 lb 1.5) ■ 10^'* fm^, the result (4J) implies a strong variation with pion mass. Such a variation 
is indeed expected from Chiral Effective Theory [|l^, |n|, |lj]. Fig. || displays the result ( pH] ) 
in relation to the pion mass dependence expected from the "Small Scale Expansion" approach 
JIO , 11], a systematic extension of leading-one-loop Heavy Baryon Chiral Perturbation Theory by 
explicit A degrees of freedom. While at a pion mass of nij^ = 759 MeV, Chiral Effective Theory 
cannot be taken as more than a qualitative hint, the comparison displayed in Fig. ||is tantalizing and 
calls for a push towards lower pion masses, where a quantitative connection with Chiral Effective 
Theory can be made. In view of the computational effort which was required for the present 
dynamical quark calculation at niji = 759 MeV, such further progress is feasible, but will require 
significant resources by current standards. 
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Figure 3: Comparison of the lattice measurement of the electric polarizability a of the neutron obtained in 
this work with the pion mass dependence expected from the "Small Scale Expansion" approach [ |l0[ [Til ]. 
This Chiral Effective Theory is expected to be accurate up to pion masses of about 300 — 400MeV (solid 
line); the dashed line is obtained if one naively continues the Chiral Effective Theory expression to higher 
pion masses, beyond its realm of controlled quantitative applicability. 
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